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ABSTRACT: We present an updated extraction of the transversity parton distribution based 
on the analysis of pion-pair production in deep-inelastic scattering off transversely polarized 
targets in collinear factorization. Data for proton and deuteron targets make it possible 
to perform a flavor separation of the valence components of the transversity distribution, 
using di-hadron fragmentation functions taken from the semi-inclusive production of two 
pion pairs in back-to-back jets in e + e _ annihilation. The e + e _ data from Belle have been 
reanalyzed using the replica method and a more realistic estimate of the uncertainties on the 
chiral-odd interference fragmentation function has been obtained. Then, the transversity 
distribution has been extracted by using the most recent and more precise COMPASS data 
for deep-inelastic scattering off proton targets. Our results represent the most accurate 
estimate of the uncertainties on the valence components of the transversity distribution 
currently available. 
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1 Introduction 

Parton distribution functions (PDFs) describe combinations of number densities of quarks 
and gluons in a fast-moving hadron. At leading twist, the spin structure of spin-half hadrons 
is specified by three PDFs. The least known one is the chiral-odd transverse polarization 
distribution h\ (transversity) because it can be measured only in processes with two hadrons 
in the initial state, or one hadron in the initial state and at least one hadron in the final 
state ( e.g . Semi-Inclusive DIS - SIDIS). 

The transversity distribution was extracted for the first time by combining data on 
polarized single-hadron SIDIS together with data on almost back-to-back emission of two 
hadrons in e + e - annihilations [1, 2]. The difficult part of this analysis lies in the fac¬ 
torization framework used to interpret the data, since it involves Transverse Momentum 
Dependent partonic functions (TMDs). QCD evolution of TMDs must be included to an¬ 
alyze SIDIS and e + e _ data obtained at very different scales, but an active debate is still 
ongoing about the implementation of these effects (see, e.g., Refs. [3-5] and references 
therein). 

Alternatively, transversity can be extracted in the standard framework of collinear fac¬ 
torization using SIDIS with two hadrons detected in the final state. In this case, h± is 
multiplied by a specific chiral-odd Di-hadron Fragmentation Function (DiFF) [6-8], which 
can be extracted from the corresponding e + e~ annihilation process leading to two back- 
to-back hadron pairs [9, 10]. In the collinear framework, evolution equations of DiFFs can 
be computed [11], Using (7T + 7T - ) SIDIS data off a transversely polarized proton target 
from HERMES [12] and Belle data for the process e + e _ —»• (7r + 7r _ )(7r + 7r _ )X [13], a point- 
by-point extraction of transversity was performed for the first time in the collinear frame¬ 
work [14], Later, including SIDIS data from transversely polarized proton and deuteron 
targets from COMPASS [15], the valence components of up and down quarks were separated 
and independently parametrized [16]. Recently, the point-by-point extraction has been ver¬ 
ified and extended also to the case of single-hadron SIDIS, showing that the transversity 
distributions obtained with two different mechanisms are compatible with each other [17]. 
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In this paper, we update the extraction of DiFFs from e + e _ annihilation data by 
performing the fit using the replica method [16]. Then, using the most recent SIDIS data for 
charged pion pairs off a transversely polarized proton target by COMPASS [18] we extract 
the transversity h\. thus obtaining the currently most realistic estimate of the uncertainties 
involved. 

In Sec. 2, we summarize the theoretical framework. In Sec. 3, we show the results of 
our updated extraction of DiFFs. In Sec. 4, we comment the salient features of the re¬ 
extracted valence components of transversity. Finally, in Sec. 5 we draw some conclusions 
and mention possible extensions of our analysis. 


2 Theoretical framework for two-hadron SIDIS 


We consider the process £(k) + N(P) —> £(k') + Hi (Pi) + ^ 2 (^ 2 ) + X, where t denotes 
the incoming lepton with four-momentum k, N the nucleon target with momentum P, 
mass M, and polarization S, Hi and H 2 the produced unpolarized hadrons with momenta 
Pi, P 2 and masses Mi, M 2 , respectively. We define the total Pf, = P\ + P 2 and relative 
R = (Pi — P2)/2 momenta of the pair, with P% = M% <C Q 2 = —q 2 > 0 and q = k — k! the 
space-like momentum transferred. As usual in SIDIS, we define also the following kinematic 
invariants 


x = 


Q 2 


2 P-q’ 


P-Ph 

Z = —— = Zi 


Z2 


y = 


C = 


P • q 
~PM ’ 
2 RP 


zi - z 2 


P-q ^ Ph-P z ’ 

where zi, z 2 , are the fractional energies carried by the two final hadrons. 


( 2 . 1 ) 

( 2 . 2 ) 



Figure 1. Kinematics of the two-hadron semi-inclusive production. The azimuthal angles <f >}j 
of the component Rt of the dihadron relative momentum , and tfis of the component St of the 
target polarization, transverse to both the virtual-photon and target-nucleon momenta q and P , 
respectively, are evaluated in the virtual-photon-nucleon center-of-momentum frame. 

The kinematics of the process is depicted in Fig. 1 (see also Refs. [12, 16]). Of particular 
relevance are the azimuthal angles of the R and S vectors. In fact, for DiFFs it is natural to 
introduce the vector Rt as the component of R perpendicular to P and P^. However, the 
cross section will depend on the azimuthal angles of both Rt and S measured in the plane 
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perpendicular to (P, q). We denote the latter ones as cpn and (f)g, respectively. In Ref. [19], 
the covariant definition of </># and 4>s is derived and compared with other non-covariant 
definitions available in the literature, pointing out the potential differences depending on 
the choice of the reference frame. For the purpose of this paper, we express <f)R and (f>s in 
the target rest frame: 


(q X k) ■ R t (qx k) • (q x R T ) 

\(q X k) ■ Rt\ \q X k\\q x Rt \ 

_ (qxk) ■ S T (q x k) ■ (q x S t ) 

| (q x k) • St\ \q x k\\q x St\ 


(2.3) 


Equation (2.3) is valid also in any frame reached from the target rest frame by a boost 
along q , up to corrections of order 0{1/Q 2 ). 

We also define the polar angle 9 which is the angle between the direction of the back-to- 
back emission in the center-of-mass (cm) frame of the two final hadrons, and the direction 
of Ph in the photon-proton cm frame (see Fig. 1). We have 

\R\ = 1 ^JmI - 2 (Mf + Mf) + (Ml - Mf )2/M h 2 , 

Rt = R sin 9 . (2.4) 


The invariant ( of Eq. (2.2) can be shown to be a linear polynomial in cos 9 [20]. 

In the one-photon exchange approximation and neglecting the lepton mass, to lead¬ 
ing order in the couplings, the differential cross section for the two-hadron SIDIS of an 
unpolarized lepton off a transversely polarized nucleon target reads [16] 

da 

dx dy ddz d<j)R dM% d cos 9 

|^(y) F uu + |Sr| B{y) sin (cj) R + <p s ) p^ ( ^ +0s) | , (2.5) 

where a is the fine structure constant, A(y) = 1 — y + y 2 / 2, B(y) = 1 — y, and the angle 
^ is the azimuthal angle of k' around the lepton beam axis with respect to the direction of 
S. In DIS kinematics, it turns out dip ~ d(f>s [21]. 

In the limit M? <C Q 2 , the structure functions in Eq. (2.5) can be written as products 
of PDFs and DiFFs [8, 20, 22]: 

Fuu = x^e 2 ff(x;Q 2 ) Dl(z, cos 9, M h \Q 2 ) , (2.6) 

<? 

p^<<t>R+<t>s) = ^2 e 2 q h\{x-, Q 2 ) H? q (z, cos 6, M h -Q 2 ) , (2.7) 

Mh q 

where e q is the fractional charge of a parton with flavor q. The Df is the DiFF describing the 
hadronization of an unpolarized parton with flavor q into an unpolarized hadron pair. The 
H^ q is a chiral-odd DiFF describing the correlation between the transverse polarization of 
the fragmenting parton with flavor q and the azimuthal orientation of the plane containing 
the momenta of the detected hadron pair. 
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Since M? <C Q 2 , the hadron pair can be assumed to be produced mainly in relative s 
or p waves, suggesting that the DiFFs can be conveniently expanded in partial waves. From 
Eq. (2.4) and from the simple relation between £ and cos0, DiFFs can be expanded in Leg¬ 
endre polynomials in cos@ [20]. After averaging over cos#, only the term corresponding to 
the unpolarized pair being created in a relative A L = 0 state survives in the D\ expansion, 
while the interference with |AL| = 1 survives for Hf [20]. The simplification holds even if 
the 6 dependence in the acceptance is not complete but symmetric about 6 = 7r/2. Without 
ambiguity, the two surviving terms will be identified with D\ and Hf, respectively. 

By inserting the structure functions of Eqs. (2.6), (2.7) into the cross section (2.5), we 
can define the single-spin asymmetry (SSA) [8, 20, 23] 


^SlDls(-T M h ;Q) 


B(y) |fl| Eq e q h-i(x; Q 2 ) Hf q (z, M h ; Q 2 ) 
Md) M h Eq egfi{x-Q 2 )Dl(z, M h -,Q 2 ) 


( 2 . 8 ) 


For the specific case of 7r + 7r _ production, isospin symmetry and charge conjugation 
suggest D q = D\ and H^ q = —H± q for q = u,d,s, and also Hf u = —H^ d and Hf 3 = 
0 [14, 16, 23]. Moreover, from Eq. (2.8) the x-dependence of transversity is more conve¬ 
niently studied by integrating the z- and M/,-dependences of DiFFs. So, in the analysis the 
actual combinations used for the proton target are [16] 


x 


h\{x-Q 2 ) = xh\ v {x]Q 2 ) - \ 

_ ^sidis(‘ t; Q") Mv) 9 
ni(Qi) B(y) 4 


xh dv (x] Q 2 ) 

e q n q(Q 2 ) X fl +q ( X :Q 2 ) , 


(2.9) 


and for the deuteron target are 

x hf(x' } Q 2 ) = x h'( v (x; Q 2 ) + xhf u (x: Q 2 ) 
_ ^sidis( x ’ Q 2 ) 


nl(Q 2 


3 \A n i(Q 2 ) + 4 n q(Q 2 )] x fi +9 ( x 'i Q 2 ) > 

q=u,d,s 


( 2 . 10 ) 


where h\ v = h\ — h\, f q+q = f q + f q . q = d,u, s if q = u, d, s, respectively ( i.e . it reflects 
isospin symmetry of strong interactions inside the deuteron), and 

n q (Q 2 ) = JdzJ dM h D q (z,M h -Q 2 ) , (2.11) 

n\(Q 2 ) = J dzj dM h ^ H? q (z, M h - Q 2 ) . (2.12) 

Using Eqs. (2.9) and (2.10), we can extract the valence components of transversity from 
the measurement of SSA Mg IDIg and Agj DIS , and from the knowledge of DiFFs through 
Eqs. (2.11) and (2.12). 


3 Extraction of Di-hadron Fragmentation Functions 

The unknown DiFFs in Eqs. (2.11) and (2.12) can be extracted from the process e"*“e - —> 
(7T + 7r - )j e t(7r + 7r - )prA. Namely, an electron and a positron annihilate producing a virtual 


- 4 - 







photon (whose time-like momentum defines the hard scale Q 2 > 0). Then, the photon 
decays in a quark and an antiquark, each one fragmenting into a residual jet. The two jets 
are produced in a back-to-back configuration; this is granted by requiring that the total 
momentum P\ of the (7r + 7r _ )jet pair in the quark jet and the total momentum P ^ of the 
(7r + 7r _ )pr pair in the antiquark jet are such that P^- P^ ~ Q 2 (in the following, all overlined 
variables will refer to the antiquark jet). 

The leading-twist cross section in collinear factorization, namely by integrating upon 
all transverse momenta but R^ and Rr, can be written as [10] 

-=—=-5?— = -—tt da 0 (1 + cos(d>R + dip) A p + p 

d cos 02dzd cos 9dMh.d(j)jidzd cos 6dMhd(j)R 47r 2 \ 



where 9-2 is the angle in the lepton plane formed by the positron direction and Py (according 
to the Trento conventions [24]), and the azimuthal angles ipn and (j)^ give the orientation of 
the planes containing the momenta of the pion pairs with respect to the lepton plane (see 
Fig.l of Ref. [10] for more details). The da 0 is the unpolarized cross section producing an 
azimuthally flat distribution of pion pairs coming from the fragmentation of unpolarized 
quarks. The term A e + e ~ represents the so-called Artru-Collins asymmetry and is given 
by [9] 


A e + e ~ ~ 


sin 2 02 
1 + COS 2 9-2 


sin# sin# 


*| R\ |R| Eg 4 H i q (z, M h - Q 2 ) H? q (z, M h -Q 2 


M h M h Eg e q D\(z, M h \ Q 2 ) D\{z, M h -, Q 2 ) 


. (3.2) 


The H^ q can be extracted from the Artru-Collins asymmetry by conveniently inte¬ 
grating upon the hemisphere of the antiquark jet. For k + tt ~) production, isospin and 
charge conjugation symmetries of DiFFs imply riq(Q 2 ) = n q (Q 2 ) for q = u,d,s,c, and 
n]j(Q 2 ) = —nl(Q 2 ) for q = u,d, neglecting other components and with the further con¬ 
straint ni,(Q 2 ) = —n^(Q 2 ). Then, Eq. (3.2) is simplified to [10] 


sin 2 #2 
1 + COS 2 02 


. Q . -5 H(z,M h ;Q 2 ) 
sin 0 sin 0 —————— 

9 D(z, M h ] Q 


2\ ’ 


where 


and 


D(z, M h - Q 2 ) = ^ D?(z, M h - Q 2 ) n u (Q 2 ) + ^ Df(z, M h ; Q 2 ) n d {Q 2 ) 
+ ^ D l(z, M h ] Q' 2 ) n s (Q 2 ) + ^ Dl(z,M h ]Q 2 )n c (Q 2 ) 


H(z, M h ] Q 2 ) = i ^H< u (z,M h ]Q 2 )nl(Q 2 ) 
My, 


1 


1 + COS 2 02 9 
sin 2 02 5 sin 9 sin 9 


D(z,M h ] Q 2 ) A e+e 


with the normalization 


j dzj dM h H(z, M h , Q 2 ) = [nl(Q 2 )] 2 


(3.3) 


(3.4) 


(3.5) 


(3.6) 
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Since a measurement of the unpolarized differential cross section is still missing, the un¬ 
polarized DiFF D\ is taken from our previous analysis in Ref. [10], where it was parametrized 
to reproduce the two-pion yield of the PYTHIA event generator tuned to the Belle kinematics. 
The fitting expression at the starting scale Qq = 1 GeV 2 was inspired by previous model 
calculations [ 8 , 23, 25, 26] and it contains three resonant channels (pion pair produced by 
p, u), and Kg decays) and a continuum. For each channel and for each flavor q = u,d, s, c, 
a grid of data in (z, Mh) was produced using PYTHIA for a total amount of approximately 
32000 bins. Each grid was separately fitted using the corresponding parametrization of D\ 
and evolving it to the Belle scale at Q 2 = 100 GeV 2 . An average x 2 P er degree of freedom 
(x 2 /d.o.f.) of 1.62 was reached using in total 79 parameters. More details can be found in 
Ref. [10] . 

As for the polarized DiFF, we deduce the experimental value of Eq. (3.5) in each bin, 
denoted H exp , by using the experimental data for the Artru-Collins asymmetry A e + e - and 
the corresponding average values of the angles 0 2 , 9, 9, taken from Ref. [13]. The function 
D is calculated from Eq. (3.4) using the unpolarized DiFFs D\ resulting from the fit of the 
PYTHIA’s two-pion yield. The experimental data for A e + e ~ are organized in a (z, Mh) grid 
of 64 bins [13]. Some of them are empty or scarcely populated; therefore, we used only 46 
of them [ 10 ]. 

The fitting value in each bin, denoted H th , is based on the following expression at the 
starting scale Qg = 1 GeV 2 [10]: 


H(z,M h ,Q 2 0 -{p}) = N h 2|i2| (1 — z) exp[ 7 l (*- 72 M ft )]BW (m p ,^-;M h ) 

V m p / 


P( 0, l,<5i, 0,0; 2 ) + zP(0, 0, 62 , S 3 , 0; M h ) + - P(0,0, £ 4 , <J 5 ,0; M h ) 


(3.7) 


where {p} denotes the vector of 9 parameters {p} = (Nh, 7 i, 7 2 , <5i, S 2 , <$ 3 , £ 4 , 65 , rj). The 
polynomial P and the Breit-Wigner function BW are defined by 

P(a 1; a 2 , 03 , 04 , 05 ; a:) = a\- + o 2 + a 3 x + a 4 x 2 + a 5 x 3 , 

x 

BW(m, T; x) = ^ 2 _ m2)2 + m 2 p 2 • (3.8) 

The function BW is proportional to the modulus squared of a relativistic Breit-Wigner for 
the considered resonant channel, and it depends on its mass m and width T. In our case, 
the p —> ( 7 r + 7 r _ ) decay involves the fixed parameters m p = 0.776 GeV and T p = 0.150 GeV. 

The function H of Eq. (3.7) is then evolved to the Belle scale using the H0PPET code [27] 
suitable extended to include chiral-odd splitting functions at leading order [10]. We have 
used two different values of a s (M^) in the evolution code, namely 0.125 [28] and 0.139 [29], 
in order to account for the theoretical uncertainties on the determination of the Aqcd 
parameter. At variance with Ref. [10], the error analysis is carried out using the same 
Monte Carlo approach adopted in our previous extraction of transversity [16]. The approach 
consists in creating N replicas of the data points. In each replica (denoted by the index 
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1.0 


R(M/i.z) 


R (M h ,z) 




Figure 2. The ratio R(z,Mh) of Eq. (3.10) as a function of M /, at Qq = 1 GeV 2 for three 
different 2 = 0.25 (shortest band), z = 0.45 (lower band at Mh ~ 1.2 GeV), and z = 0.65 (upper 
band at Mh ~ 1.2 GeV). Left panel for results obtained with a s (M§) = 0.125, right panel with 
a s (M§) = 0.139. For the calculation of the uncertainty bands see details in the text. 


r), the data point in the bin ( Zi,Mhj ) is perturbated by a Gaussian noise with the same 
variance as the experimental measurement. Each replica, therefore, represents a possible 
outcome of an independent measurement; for the bin ( Zi,Mhj ), we denote it by 
The number of replicas is chosen N = 100 in order to accurately reproduce the mean and 
standard deviation of the original data points. The standard minimization procedure is 
applied to each replica r separately, by minimizing the following error function 


*?({p» = E 

ij 



(3.9) 


where H X A is the fitting value of Eq. (3.7) depending on the vector {p} of 9 parameters, 
and the error (jjj for each replica is taken to be equal to the error on the original data point 
for the bin (, Zi,Mhj ). In fact, in the expression of H from Eq. (3.5) the dominant source 
of uncertainty comes from the experimental error on the measurement of A e + e -. The very 
large statistics available from PYTHIA in the Monte Carlo simulation of two-pion yields makes 
the statistical uncertainty on D negligible [10]. However, there is still a source of systematic 
error that is not taken into account in this analysis. This problem can be overcome only 
when real data for the unpolarized cross section will become available. Meanwhile, in this 
work <7 l j is obtained by summing in quadrature the statistical and systematic errors for the 
measurement of A e + e - reported by the Belle collaboration [13], multiplied by all factors 
relating A e + e - to H according to Eq. (3.5). 

The vector of parameters that initializes the minimization in Eq. (3.9) corresponds to 
the one that produced the best fit in the previous analysis of Ref. [10] using the standard 
Hessian method. Then, the minimization results in N different vectors of best-fit param¬ 
eters, {po}> v = 1 These vectors can be used to produce N different values for 

any theoretical quantity. Using Eqs. (3.5) and (3.6), we can produce N different replicas 
of the polarized DiFF Hf u . The N theoretical outcomes can have any distribution, not 
necessarily Gaussian. Hence, the la confidence interval is in general different from the 68 % 
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1.0 


R(M/i.z) 


R (M h ,z) 




Figure 3. The ratio R(z,Mh) of Eq. (3.10) as a function of z at Qq = 1 GeV 2 for three different 
Mh = 0.4 GeV (lower band at z ~ 0.8), Mh = 0.8 GeV (mid band at z ~ 0.8), and Mh = 1 GeV 
(upper band at z ~ 0.8). Left panel for results obtained with a fl (M§) = 0.125, right panel with 
a s (M |) = 0.139. For the calculation of the uncertainty bands see details in the text. 


interval which, in our case, can simply be obtained by rejecting for each experimental point 
( Zi,Mhj ) the largest and the lowest 16% of the N values. This approach produces a more 
realistic estimate of the statistical uncertainty on DiFFs. In fact, we noticed that the mini¬ 
mization often pushes the theoretical functions towards their upper or lower bounds, where 
the x 2 does no longer display a quadratic dependence upon the parameters. Instead, the 
Monte Carlo approach does not rely on the prerequisites for the standard Hessian method 
to be valid. Although the minimization is performed on the function defined in Eq. (3.9), 
the agreement of the N theoretical outcomes with the original Belle data is better expressed 
in terms of the standard % 2 function [30]. The x 2 can be obtained by replacing H^ v r in 
Eq. (3.9) with the corresponding value inferred from the original data set without Gaussian 
noise. 

We show our results through the following ratio: 


R(z, M h ) 


|fl| H< u {z,M h -Ql) 
M h Df(z,M h -Ql) ’ 


(3.10) 


where both DiFFs are summed over all fragmentation channels and the ratio is evaluated 
at the hadronic scale Q'f } = I GeV 2 . In Fig. 2, we consider the ratio R as a function 
of the invariant mass Mh for three different values of the fractional energy z: z = 0.25 
(shortest band), 2 = 0.45 (lower band at Mh ~ 1.2 GeV), and z = 0.65 (upper band at 
M h 1.2 GeV). The left panel displays the results with a s (M|) = 0.125, the right one 
with a S (M|) = 0.139. Each band corresponds to the 68% of all N = 100 replicas, produced 
by rejecting the largest and lowest 16% of the replicas’ values for each ( z,Mh ) point. The 
shortest band (for 2 = 0.25) stops around Mh ~ 0.9 GeV because there are no experimental 
data at higher invariant masses for such low values of z. In this kinematic range, the fit is 
much less constrained and, consequently, the uncertainty band becomes larger. Comparing 
the two panels reveals a mild sensitivity of the results to the choice of a s (M§), hence of 
Aqcd- Figure 2 represents an update of the upper panel of Fig. 6 in Ref. [10], with a more 
realistic estimate of the statistical uncertainty on the polarized DiFF H^. 
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Figure 4. Histogram of the distribution of N = 100 x 2 /d.o.f. for a s [M\) = 0.125 (left panel) and 
a s (M§) = 0.139 (right panel). The solid curve corresponds to a Gaussian distribution centered 
on the average of the N = 100 x 2 values. The shaded area represents the ler variance. The 
normalization of the Gaussian distribution is adapted to the histogram profile. 


In Fig. 3, the same quantity R of Eq. (3.10) is plotted as a function of z for three 
different values of the invariant mass: M/j = 0.4 GeV (lower band at z ~ 0.8), Mh = 0.8 
GeV (mid band at z ~ 0.8), and = 1 GeV (upper band at z ~ 0.8). Again, the left 
panel displays the results with a s (M|) = 0.125 while the right one with a s (M |) = 0.139. 
The results of the two panels are similar, as we found in the previous figure, except for the 
bands with larger M\ at very high values of 2 . Again, Fig. 3 represents a realistic update 
of the bottom panel of Fig. 6 in Ref. [10]. 

From the minimization in Eq. (3.9), we obtain N different x 2 > each one corresponding 
to a different vector of fitting parameters {p} r , r = 1... N. If the model is able to give a 
good description of the data, the distribution of the N values of x 2 /d.o.f. should be peaked 
at around one. In real situations, the rigidity of the model shifts the position of the peak to 
higher values. In Fig. 4, we show the histogram for the distribution of the N values of the 
X 2 /d.o.f. It is not peaked at 1 but slightly above 1.5. For sake of illustration, we compare it 
with the solid line representing a Gaussian distribution centered at the average value of the 
N different x 2 /d.o.f. The shaded area represents its la variance. As in previous figures, 
the left panel shows the x 2 /d.o.f. obtained using a s (M|) = 0.125 in the evolution code, 
while in the right panel a s (M§) = 0.139. The salient features of the x 2 distribution remain 
the same in both cases. 

Finally, in Tab. 1 we show the average value {(p)} and variance {a p } for each of the 
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a s (Mf) = 0.125 

a a (Af§) = 0.139 

M 

{<P» 

Wp} 

«P>} 

{o-p} 

N 

0.016 

0.008 

0.015 

0.007 

7i 

-2.73 

0.91 

-2.27 

0.85 

72 

-0.92 

0.84 

-1.10 

0.85 

O 

23 

21 

21 

23 

<52 

-200 

48 

-195 

53 

<53 

278 

30 

268 

34 

<54 

39 

11 

43 

13 

<5s 

-44 

13 

-48 

13 

V 

0.29 

0.12 

0.30 

0.09 


Table 1 . Average value {(p)} and variance {cr p } of each element of the vector of fitting parameters 
{p} in Eq. (3.7), calculated by fitting N = 100 replicas of the experimental data points for the Artru- 
Collins asymmetry (see discussion in the text around Eq. (3.9)). Left columns: QCD evolution 
performed with a s (M§) = 0.125; right columns with a s (M = 0.139. 


9 elements in the vector {p} of the fitting parameters in Eq. (3.7). They are calculated 
from the set of N different values obtained by fitting the N replicas of the experimental 
data points, i.e. by minimizing the function E P ({p}) in Eq. (3.9) for r = The 

left pair of columns of numbers are obtained using a s {M'%) = 0.125, the right ones using 
a s (Mz) = 0.139. For some parameters, the variance is large compared to their average 
value, indicating that they are loosely constrained by the fit although the resulting \ ,2 /d.o.f. 
are reasonable (see Fig. 4). 

4 Extraction of transversity 

The valence components of transversity are extracted by combining Eqs. (2.9) and (2.10). 
In these equations, there are three main external ingredients: the unpolarized distributions 
ff, the single-spin asymmetries ^TglDlS ’ an< ^ the integrals n q and n£. The unpolarized 
distributions /[ ; are taken from the MSTW08 set [29] at leading order (LO). 

The experimental data for the single-spin asymmetries Ag IDIg and AgJ DIg are taken 
from the HERMES and COMPASS measurements on di-hadron SIDIS production off trans¬ 
versely polarized proton and deuteron targets. In our previous extraction [16], we used the 
HERMES data for a proton target from Ref. [12], and the COMPASS data for unidentified 
hadron pairs h + h~ produced off deuteron and proton targets from Ref. [15] (corresponding 
to the 2004 and 2007 runs, respectively). In an intermediate step [31], we have updated our 
extraction by using the new COMPASS data for unidentified hadron pairs h + h~ produced 
off protons [32], corresponding to the 2010 run. In this work, we select the new COMPASS 
data for identified vr + 7r~ pairs produced off proton targets [18], still corresponding to the 
2010 run. More precisely, we have used the results of Tab. A. 10 and A.26 presented in the 
appendix A.4 of Ref. [33] where the proton data of the 2010 run have been complemented 
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with a homogeneous re-analysis of all identified pairs from the 2004 and 2007 runs using 
the same data selection, methods, and binning. 

Finally, the third ingredient is represented by the integrals n q and nt of Eqs. (2.11) 
and (2.12), respectively. They are evaluated according to the appropriate experimental 
cuts: 0.2 < z < 1 and 0.5 GeV < M\ < 1 GeV for HERMES, 0.2 < z < 1 and 0.29 GeV < 
Mh < 1.29 GeV for COMPASS. The arguments of the integrals, namely the DiFFs Hf u 
and D\ with q = u,d, s, c, are determined along the lines described in the previous section. 

Equations (2.9) and (2.10) are fitted using the same strategy adopted in our previ¬ 
ous extraction of Ref. [16], and here used also for the extraction of DiFFs (see previous 
section). Namely, for each bin ( Xi,Qf) a set of M replicas of the data points A giDig and 
A^dis i s created introducing a Gaussian noise with the same variance as the correspond¬ 
ing experimental measurement. Using the above information for the other ingredients ff, 
n q and rt) u . we build M different replicas of Eqs. (2.9) and (2.10) that we indicate with 
x i hi rexp (xi, Qf) and X{ h^ rexp (xj, Qf), respectively, with r = 1,... M. Again, we checked 
that with M = N = 100 replicas the mean and standard deviation of the original data 
points are accurately reproduced. In principle, for each bin (xi,Qf) we can arbitrarily se¬ 
lect one out of the N possible values of n q (Qf) and nl(Q 2 ) obtained from the N different 
replicas of DiFFs at that scale. So, in order to build Xihff/I? (xi,Qf) we simply pick up 
the replica r of DiFFs, we calculate the corresponding integrals at the scale Qf, and we 
associate them to the replica r of the data points A g jp ig in the bin (xj, Qf)- Then, for each 
replica r we separately minimize an error function similar to Eq. (3.9): 


K 2 ({p}) = E 


[®i Qb {P'})~ x i K'fexp (®i. Qi) 


Ah 


p/D 
1 data 


(xi,Qi)) 


(4.1) 


where Ahf / d f ta ( x i,Qf) are the errors on the original data points Agj^ ig multiplied by all 
factors according to Eqs. (2.9) and (2.10), in close analogy with aij in Eq. (3.9). 

The function Xi (xi, Qf, {p'}) in Eq. (4.1) is obtained by evolving at the scale Qf 

of each bin a fitting function that depends on the vector of parameters {p'}- The main 
theoretical constraint that transversity must satisfy at each scale is Soffer’s inequality [34], 
If it is verified at some initial Qf, it will do also at higher Q 2 > Qf [35]. Therefore, following 
our previous work [16] we have parametrized the fitting function for the valence flavors u v 
and d v at Qf = 1 GeV 2 as 


x h\ v (x; Qf) = tanh x 1//2 (A q + B q x + C q x 2 + D q x 3 ) x SB 9 (x; Qq) + SB 9 (x; Qq) 


(4.2) 


where the analytic expression of the Soffer bound SB' ? (x; Q 2 ) can be found in the Appendix 
of Ref. [16]. The implementation of the Soffer bound depends on the choice of the unpo¬ 
larized PDF, as mentioned above, and of the helicity PDF that we take from the DSSV 
parameterization [36]. As in Ref. [16], the error on the Soffer bound is dominated by the 
uncertainty on the helicity <j\, that we checked to be negligible with respect to the experi¬ 
mental errors on Agj^ ig . In Eq. (4.2), the hyperbolic tangent is such that the Soffer bound 
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x h“'(x)+x hj v (x) 



Figure 5. The combinations of Eq. (2.9), left panel, and Eq. (2.10), right panel. The black circles 
are obtained from the HERMES data for the SSA Ag IDIg ; the lighter squares from the COMPASS 
data for both Ag IDIS and A^? IDI g. The uncertainty band represents the selected 68 % of all fitting 
replicas in the rigid scenario with a s (M§) = 0.125 (see text). 


is always fulfilled. The low-x behavior is determined by the x 1 / 2 term, which is imposed 
by hand to grant that the resulting tensor charge is finite. Present fixed-target data do not 
allow to constrain it. This choice has little effect on the region where data exist, but has a 
crucial influence on extrapolations at low x. The functional form in Eq. (4.2) is very flexible 
and can display up to three possible nodes. In analogy with Ref. [16], we have explored 
three different scenarios: a) the rigid scenario, where C q = D q = 0 and the vector of param¬ 
eters {p'} contains only 4 free parameters; b) the flexible scenario, with D q = 0 and 6 free 
parameters; c) the extraflexible scenario, with all possible 8 free parameters. The function 
x* tti Qi > {'P'Y) i n Eq. (4.1) is then built by taking the proper flavor combinations of 
the fitting function in Eq. (4.2) and evolving them to the scale Qf of each bin i. Evolution is 
realized using the HOPPET code, suitably extended to include chiral-odd splitting functions 
at leading order [10], and with both values of a s (M^) = 0.125 and a s (M ]|) = 0.139. 

In Fig. 5, the points represent the combinations of Eq. (2.9) in the left panel, and of 
Eq. (2.10) in the right panel. The error bars are mainly determined by the experimental 
errors on M gIDIg and Ag^ DIg , respectively, because the uncertainty on the extracted DiFFs 
is much smaller. The black circles in the left panel are obtained when using for M gIDIg the 
HERMES measurement from Ref. [12]. The lighter squares in both panels correspond to 
the COMPASS measurements of M gIDIg (left) and Mg[ DIg (right) from Ref. [33], as explained 
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Figure 6. Histogram of the distribution of M = 100 x 2 /d.o.f. when minimizing Eq. (4.1) for 
a s (M%) = 0.125 and in the rigid scenario. The solid curve corresponds to a Gaussian distribution 
centered at the average of the M = 100 x 2 values. The shaded area represents the lcr variance. 
The normalization of the Gaussian distribution is adapted to the histogram profile. 


above. The uncertainty bands show the result of the 68% of all fitting replicas in the rigid 
scenario with a s {M z ) = 0.125. They are obtained by minimizing the error function in 
Eq. (4.1) and by further rejecting the largest 16% and the lowest 16% of the M = 100 
replicas’ values in each x point. 

In Fig. 6, we show the histogram for the distribution of the M values of the x 2 /d.o.f 
obtained by minimizing the error function in Eq. (4.1) for the rigid scenario with a s {M |) = 
0.125. For sake of illustration, we compare it with the solid line representing a Gaussian 
distribution centered around the average 1.42 of the x 2 /d.o.f. values for this scenario. The 
shaded area represents the lcr variance. The distribution is not peaked at 1 but around 
1.4 because of the rigidity of the fitting model. When changing evolution parameter from 
a s (Mg) = 0.125 to a s (M z ) = 0.139, the salient features of the x 2 distribution remain 
substantially the same and the average x 2 /d.o.f. increases by less than 3%, as it can be 
realized by inspecting Tab. 2. 


X 2 / d.o.f. 

a a (M§) = 0.125 

a s {M 2 z ) = 0.139 

rigid 

1.42 

1.46 

flexible 

1.65 

1.71 

extraflexible 

1.97 

2.07 


Table 2. The average x 2 /d.o.f. obtained by minimizing the error function in Eq. (4.1) for the three 
different scenarios explored in the fitting function, and for the two values of a s in the evolution 
code. 


In Fig. 7, we show the up valence transversity, as a function of x at Q 2 = 2.4 

GeV 2 in the flexible scenario. The brightest band in the background with dashed borders 
is the 68% of all replicas from our previous extraction [16]. The light grey band in the 
foreground with dot-dashed borders shows the 68% of all replicas obtained in this work 
when using a s (M ]|) = 0.139. The darkest band with solid borders is the result when 
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Figure 7. The up valence transversity as a function of a; at Q 2 = 2.4 GeV 2 in the flexible scenario. 
The brightest band in the background with dashed borders is the 68% of all replicas from our 
previous extraction [16]. The light grey band in the foreground with dot-dashed borders is the 68% 
of all replicas obtained in this work with a s (M§) = 0.139. The darkest band with solid borders is 
the same but for a s (M§) = 0.125. The thick solid lines indicate the Soffer bound. 


using a s (Mz) = 0.125. Finally, the thick solid lines indicate the Soffer bound. The fact 
that the latter two bands overlap almost completely confirms that our new extraction is 
not very sensitive to the value of a s (M'%), namely to the theoretical uncertainty in the 
evolution equations. On the other side, the impact of the new COMPASS data is rather 
evident. There is still overlap between present and previous extractions, but the better 
statistical precision of data produces a narrower uncertainty band, at least in the range 
0.0065 < x < 0.29 where there are data. Moreover, the replicas spread out over values that 
on average are smaller than before. Since the new COMPASS analysis of Ref. [32] deals with 
proton targets, the combination in Eq. (2.10) is not affected. Our extraction of the down 
valence transversity is basically unchanged with respect to the previous one [16]; therefore, 
we will not show it. Similar results are obtained when switching to other scenarios in the 
fitting function; we will not show them as well. 

In Fig. 8 , we show how our new results compare with other extractions of transversity 
based on the Collins effect. In the left (right) panel, the up (down) valence transversity is 
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Figure 8 . The up (left) and down (right) valence transversities as functions of a; at Q 2 = 2.4 
GeV 2 . The darker band with solid borders in the foreground is our result in the flexible scenario 
with a s (M§) = 0.125. The lighter band with dot-dashed borders in the background is the most 
recent transversity extraction from the Collins effect [2]. The central thick dashed line is the result 
of Ref. [5]. The thick solid lines indicate the Softer bound. 

displayed as a function of x at Q 2 = 2.4 GeV 2 . The darker band with solid borders in the 
foreground is our result in the flexible scenario with a s (M'%) = 0.125. The lighter band 
with dot-dashed borders in the background is the most recent transversity extraction of 
Ref. [2] using the Collins effect but applying the standard DGLAP evolution equations only 
to the collinear part of the fitting function. The central thick dashed line is the result of 
Ref. [5], where evolution equations have been computed in the TMD framework. 

In the right panel, the disagreement between our result for xhf"(x) at x > 0.1 and 
the outcome of the Collins effect is confirmed with respect to our previous analysis (see 
Fig. 4 in Ref. [16]). This is due to the fact that the COMPASS data for Ag^ DIg off deuteron 
targets remain the same. This trend is confirmed also in the other scenarios, indicating 
that it is not an artifact of the chosen functional form. As a matter of fact, our replicas for 
the valence down transversity tend to saturate the lower limit of the Soffer bound because 
they are driven by the COMPASS deuteron data, in particular by the bins number 7 and 
8 . It is worth mentioning that some of the replicas outside the 68% band do not follow 
this trend. Their trajectories are spread over the whole available space between the upper 
and lower limits of the Soffer bound, still maintaining a good x 2 /d.o.f. (typically, around 
2). It is also interesting to remark that the dashed line from Ref. [5], although in general 
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agreement with the other extraction based on the Collins effect, also tends to saturate the 
Soffer bound at x > 0.2. 

Apart from the range x > 0.1, there is a general consistency among the various extrac¬ 
tions which is confirmed also for the valence up transversity (left panel), at least for the 
range 0.0065 < x < 0.29 where there are data. This is encouraging: while the dihadron 
SIDIS data are a subset of the single-hadron ones, the theoretical frameworks used to in¬ 
terpret them are very different. Nevertheless, we point out that the collinear framework, in 
which our results are produced, represents a well established and robust theoretical context. 
On the contrary, the implementation of the QCD evolution equations of TMDs needed in 
the study of the Collins effect still contains elements of arbitrariness (see Refs. [3-5] and ref¬ 
erences therein). Moreover, we believe that our error analysis, based on the replica method 
applied to the extraction of both the DiFFs from e + e _ data and the transversity from SIDIS 
data, represents the current most realistic estimate of the uncertainties on transversity. It 
also clearly shows that we have no clue on the transversity for large x > 0.3 where there are 
no data at present. This is particularly evident in the left panel of Fig. 8: the replicas in 
the darker band tend to fill all the available phase space within the solid lines of the Soffer 
bound, graphically visualizing our poor knowledge of xh™ v (x) in that range. Similarly, data 
are missing also for very small x. and this prevents from fixing the behaviour of transversity 
for x — y 0 in a less arbitrary way than the choice made in Eq. (4.2). 

In Fig. 9, we show the "truncated" tensor charge 

/■^max 

5^q v (Q 2 ) = dxhf(x,Q 2 ), (4.3) 

^ x min 

namely the truncated first Mellin moment of the valence transversity. The integral is 
computed for x m i n = 0.0065 < x < x max = 0.29, i.e. in the range of experimental data, 
thus avoiding any numerical uncertainty produced by extrapolation outside this range. In 
the left panel, we show 5^u v (Q 2 = 10 GeV 2 ), in the right panel 6 ^d v (Q 2 = 10 GeV 2 ). They 
are calculated at Q 2 = 10 GeV 2 in order to compare with the results of Ref. [5], which 
are indicated in both panels by the leftmost circle with label 2. The black squares with 
labels 3-5 indicate our result with a s (M ]|) = 0.125 for the rigid, flexible, and extraflexible 
scenarios, from left to right respectively. The triangles with labels 6-8 correspond to the 
choice a s (Mz) = 0.139 in the same order. The corresponding error bars are computed by 
considering the distance between the minimum and the maximum values of the 68% of all 
replicas; the squares and triangles identify their equidistant point. Our results are basically 
insensitive to the choice of a s \ so, in the following we will show results only for the choice 
a s (Mg) = 0.125, forwarding the reader to Tab. 3 for the numerical values of all considered 
cases. 

In Fig. 10, we show the full Mellin moments of valence transversity at Q‘q = 1 GeV 2 , 
i.e. the tensor charges 

dq v {Q 2 ) = [ dx h\ v (x, Q 2 ) . (4.4) 

Jo 

The integration is now extended to the full x domain by extrapolating h^ v (x) outside the 
experimental range. As in the previous figure, the left panel refers to the valence up quark 
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Figure 9. Truncated tensor charges (see text) at Q 2 = 10 GeV 2 for the valence up (left panel) 
and down quark (right panel). From left to right: circle (label 2) for the value obtained through 
the Collins effect in Ref. [5], black squares (labels 3-5) for the rigid, flexible, extraflexible scenarios, 
respectively, here explored with a s (M§) = 0.125, triangles (labels 6-8) for the corresponding ones 
with a s (Mz) = 0.139. 


while the right one to the valence down quark. The two leftmost circles (labels 1, 2) are the 
results obtained from the analysis of the Collins effect using two different methods for the 
extraction of the Collins function from e + e~ annihilation data [2], The three black squares 
(labels 3-5) correspond to the results of our previous analysis [16] for the rigid, flexible, and 
extraflexible scenarios, from left to right respectively. The three rightmost triangles (labels 
6 - 8 ) indicate the outcome of the present work with a s (M§) = 0.125 in the same order. 
Consistently with Fig. 7, our new results for the up quark are smaller than the previous 
ones. They also appear globally in better agreement with the values from Ref. [2] (and 
not far from the ones obtained from the parametrization of chiral-odd Generalized Parton 
Distributions of Ref. [37]), although the large uncertainties introduced by the numerical 
extrapolation smooth most of the differences. This is particularly evident for the down 
quark, where in addition the numerical values are very close because the experimental data 
for Tigris are the same as before. 

In Fig. 11, we show the isovector nucleon tensor charge qt = 5u v —5d v . While there is no 
elementary tensor current at tree level in the Standard Model, the nucleon matrix element 
of the tensor operator can still be defined (for a review, see Ref. [38] and references therein). 
The gx belongs to the group of isovector nucleon charges that are related to flavour-changing 
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Figure 10. Tensor charges at Qq = 1 GeV 2 for the valence up (left panel) and down quark (right 
panel). From left to right: circles (label 1, 2) for the values obtained through the Collins effect 
in Ref. [2], black squares (labels 3-5) for the rigid, flexible, extraflexible scenarios explored in our 
previous extraction of Ref. [16], triangles (labels 6 - 8 ) for the present work with a s (M\) = 0.125. 


processes. A determination of these couplings may shed light on the search of new physics 
mechanisms that may depend on them [39-42] , or on direct dark matter searches [43] . The 
vector charge gy , axial charge qa , and induced tensor charge gp. are fixed by baryon number 
conservation, neutron /3-decay, and nucleon magnetic moments, respectively [44], Also the 
pseudoscalar charge gp is, to some extext, constrained by low-energy nir + scattering [45]. 
The other isovector nucleon couplings, including gp, have been determined so far only with 
lattice QCD. 

In Fig. 11, the leftmost light square with label 1 is our new result for gp = 0.81 + 0.44 
at Q 2 = 4 GeV 2 for the flexible scenario with a s (M= 0.125. We compare it with various 
lattice computations. From left to right, the black square refers to the lattice simulation of 
RQCD at m n ~ 150 MeV with n/ = 2 NPI Wilson-clover fermions [46], the black triangle 
to that of RBC-UKQCD at m n = 330 MeV with n/ = 2 + 1 domain wall fermions [47], the 
black circle to that of LHPC at m n ~ 149 MeV with n/ = 2 + 1 HEX-smeared Wilson- 
clover fermions [48], the black inverted triangle to that of PNDME at m n = 220 MeV with 
Wilson-clover fermions on a HISQ staggered nj = 2 + 1 + 1 sea [49], the black diamond 
and star to that of ETMC at physical m n with rif = 2 twisted mass fermions and at 
m n = 213 MeV with rif = 2 + 1 + 1 twisted mass fermions, respectively [50]. Our result is 
obviously compatible with the various lattice simulations because of the very large error. 
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Figure 11. Isovector tensor charge 6 u v — Sd v at Q 2 = 4 GeV 2 . From left to right: light square 
(label 1) is our result for the flexible scenario with a s (M\) = 0.125; black square for the lattice 
result of Ref. [46] (RQCD); black triangle from Ref. [47] (RBC-UKQCD); black circle from Ref. [48] 
(LHPC); black inverted triangle from Ref. [49] (PNDME); black diamond and star from Ref. [50] 
(ETMC) with 2+1 and 2+1+1 flavors, respectively. 


As already remarked, this originates from the fact that the integral in Eq. (4.4) involves 
the extrapolation of transversity outside the x range of experimental data. From Fig. 7 it is 
evident that the replicas tend to take all values within the Soffer bounds for x > 0.3 where 
there are no data, thus increasing the uncertainty. Moreover, we stress again that there is 
also a source of systematic error related to the power x 1//2 in the fitting form of Eq. (4.2). 
The absence of data at very low x leaves this choice basically unconstrained, whereas the 
value of the integral in Eq. (4.4) heavily depends on it. 

Finally, in Tab. 3 we collect all numerical values that we have obtained for the (trun¬ 
cated) tensor charge. In the upper part of the table, we show the truncated tensor charge 
5 1 q v of Eq. (4.3) at Q 2 = 10 GeV 2 for valence up and down quarks in the rigid, flexi¬ 
ble, extraflexible scenarios for the fitting function of Eq. (4.2) with a s (M j|) = 0.125 or 
a s (M|) = 0.139 in the evolution code. In the lower part of the table, we show the results 
for the same cases but for the tensor charge 5q v of Eq. (4.4) at the starting scale Q 0 = 1 
GeV 2 . 
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S"q v (Q 2 = 10 GeV 2 ) 

valence up 

valence down 


a s {M 2 ) = 0.125 

a s (M 2 ) = 0.139 

a s (M§) = 0.125 

a s (M 2 ) = 0.139 

rigid 

0.26 ±0.05 

0.24 ±0.05 

-0.19 ±0.10 

-0.19 ±0.10 

flexible 

0.25 ±0.05 

0.24 ±0.04 

-0.25 ±0.12 

-0.24 ±0.11 

extraflexible 

0.27 ±0.05 

0.25 ±0.06 

—0.24 ± 0.11 

-0.22 ±0.10 


5q v (Ql = 1 GeV 2 ) 




rigid 

0.49 ± 0.09 

0.43 ± 0.08 

0.05 ±0.25 

0.04 ±0.24 

flexible 

0.39 ±0.15 

0.40 ±0.14 

-0.41 ±0.52 

-0.32 ±0.51 

extraflexible 

0.35 ±0.14 

0.36 ±0.12 

-0.04 ±0.77 

-0.12 ±0.74 


Table 3. Summary of numerical values for the tensor charge. Upper part for the truncated 
tensor charge of Eq. (4.3) at Q 2 = 10 GeV 2 for valence up and down quarks in the rigid, flexible, 
extraflexible scenarios for the fitting function of Eq. (4.2) with a s (M|) = 0.125 or a s (M%) = 0.139 
in the evolution code. Lower part for the tensor charge of Eq. (4.4) at Qq = 1 GeV 2 . 

5 Conclusions 

The transversity parton distribution function is an essential piece of information on the 
nucleon at leading twist. Its first Mellin moment is related to the nucleon tensor charge. 
Due to its chiral-odd nature, transversity cannot be accessed in fully inclusive deep-inelastic 
scattering (DIS). Within the framework of collinear factorization, it is however possible 
to access it in two-particle-inclusive DIS in combination with Dihadron Fragmentation 
Functions (DiFFs). The latter can be extracted from e + e~ annihilations producing two 
back-to-back hadron pairs, and evolution equations are known to connect DiFFs at the 
different scales of the two reactions. 

In this paper, we have updated our first extraction of DiFFs from e + e _ annihilation 
data in Ref. [10] by performing the error analysis with the so-called replica method. The 
method is based on the random generation of a large number of replicas of the experimental 
points, in this case the Belle data for the process e + e _ —> (7r + 7r - ) (7r + 7r~) [13]. Each replica 
is then separately fitted, producing an envelope of curves whose width is the generalization 
of the lcr uncertainty band when the distribution is not necessarily a Gaussian. As such, 
this method allows for a more realistic estimate of the uncertainty on DiFFs. 

As a second step, we have used the above result to update our first extraction of the 
up and down valence transversities in a collinear framework [16], employing data for two- 
particle-inclusive DIS off transversely polarized proton and deuteron targets. In particular, 
we have considered the recent measurement from the COMPASS collaboration for identified 
hadron pairs produced off transversely polarized proton targets [18]. We have randomly 
generated replicas of these data and we have fitted them, making an error analysis similar 
to what has been done for DiFFs. We have noticed that many of these trajectories hit the 
Soffer bound, i.e. they lie close to the borders of the phase space where the \' 2 function 
cannot be expected to have the quadratic dependence on the fit parameters as required by 
the standard Hessian method. Hence, we stress again that the replica method allows for a 
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more reliable error analysis and we believe that the results shown in this paper represent 
the currently most realistic estimate of the uncertainties on the transversity distribution. 


As in our previous extraction [16], we have adopted different scenarios for the functional 
form, all subject to the Soffer bound. We have further explored the sensitivity to the 
theoretical uncertainty on Aqcd by using two different prescriptions for a s (Af|) [28, 29]. 
In the range of experimental data, our results show little sensitivity to the variation of these 
parameters. The 68% band for the valence up transversity turns out to be narrower than in 
the previous extraction because of the more precise COMPASS data on the proton target. 
Nevertheless, there is a significant overlap with the other existing parametrizations based 
on the Collins effect in single-hadron-inclusive DIS [2, 5]. The only source of discrepancy 
lies in the range x > 0.1 for the valence down quark, where all replicas are driven to hit the 
lower Soffer bound irrespectively of the functional form and evolution parameters adopted. 
This behavior is induced by two specific bins in the set of COMPASS experimental data 
for the deuteron target. Since this data set is not changed with respect to our previous 
extraction, the present results just confirm those findings in Ref. [16]. It is interesting to 
note that also the down transversity of Ref. [5], extracted from the Collins mechanism but 
with evolution effects described in the TMD framework, tends to saturate the Soffer bound 
at x > 0.2. 


We have also calculated the first Mellin moment of transversity, i.e. the tensor charge, 
either by computing the integral upon the range of experimental data (truncated tensor 
charge) or by extrapolating the transversity to the full support [0,1] in the parton fractional 
momentum x. The latter option obviously induces a much larger error that somewhat 
decreases the relevance of the observed compatibility with the results obtained from the 
extraction based on the Collins effect. Nevertheless, we find good agreement also for the 
truncated tensor charges obtained in Ref. [5]. We have also computed the isovector tensor 
charge qt- The determination of the latter may shed light on hypothetical new elementary 
electroweak currents that are being explored through neutron f3 decays [39-42], or even 
on direct dark matter searches [43]. Our result has a very large error because, again, it 
requires the extrapolation of transversity outside the range of experimental data. Anyway, 
it is compatible with all the lattice results available in the literature. 


The large uncertainty caused by extrapolating the transversity reflects the need of 
two-particle-inclusive DIS data either at large and at very small x. More data have been 
released by the COMPASS collaboration that include also different types of hadron pairs 
(e.g., Kit) [18] and should allow to improve the flavour separation of transversity. More 
insight along the same direction will come also from polarized proton-proton collisions [51], 
where data for the semi-inclusive production of hadron pairs are expected from the PHENIX 
and STAR collaborations (see, e.g., Ref. [52]). Finally, two-particle-inclusive DIS will be 
measured also at JLab in the near future, which should considerably increase our knowledge 
of transversity at large x. 
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